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receive input corresponding to user-selected ^\ q 
point near nodule boundary 



pre-process region-of-interest around point 
using Expectation Maximization 
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construct model of nodule from point, 
to initalize a set of parameters \l that 230 
represent nodule shape 
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obtain estimates of continuous parts and 
discontinuities of boundary, and update 
estimates until (I converges, using dynamic 240 
programming and Expectation Maximization 
(EM) 
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segment nodule, based on estimates of 
continuous parts and discontinuities of 250 
boundary 
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receive input corresponding to user-selected 
point near nodule boundary 
I 

pre-process region-of-interest around point 
using Expectation Maximization 

construct model of nodule from point, 
to initalize a set of parameters R>, 
where ft = [O, s]T 

I 

represent realization B of nodule boundary 
by a sum B= (U, 5 C! ) U(U,- , where 
B ci represents continuous parts and 
B.. represents discontinuities of boundary 



estimate B . and B 
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using ft, divide boundary into B ci 
and 5^. , according to criterion based 
on image data 
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estimate B ci , using MAP estimate 
according to 

B ci = arg maxB ci p(B ci \I, E) _ 

/ being a slice from 3D CT dataset 
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estimate MAP density as 
p(Bcixl, R>) = 1/z exp(E&(Bci)), where 
Eft(Bci) is the sum of internal shape 
and external image energies 
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minimize Efi(Bci), using time-delayed 
discrete dynamic programming 
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connect estimated B ri , to obtain an 



estimate on B 
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update ft, using circle fitting method 
and estimated B. and B A 



segment nodule, based on estimated 

B ci and B dj 
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